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Abstract. 

Accurate investigations of quantum level energies in molecular systems are shown 
to provide a test ground to constrain the size of compactified extra dimensions. 
This is made possible by the recent progress in precision metrology with ultrastable 
lasers on energy levels in neutral molecular hydrogen (H 2 , HD and D 2 ) and the 
molecular hydrogen ions (H^, HD+ and D^). Gomparisons between experiment and 
quantum electrodynamics calculations for these molecular systems can be interpreted 
in terms of probing large extra dimensions, under which conditions gravity will become 
much stronger. Molecules are a probe of space-time geometry at typical distances 
where chemical bonds are effective, i.e. at length scales of an A. Constraints on 
compactification radii for extra dimensions are derived within the Arkani-Hamed- 
Dimopoulos-Dvali framework, while constraints for curvature or brane separation are 
derived within the Randall-Sundrum framework. Based on the molecular spectroscopy 
of D 2 molecules and HD"'' ions, the compactification size for seven extra dimensions (in 
connection to M-theory defined in 11 dimensions) of equal size is shown to be limited 
to R^ < 0.6/rm. While limits on compactification sizes of extra dimensions based on 
other branches of physics are compared, the prospect of further tightening constraints 
from the molecular method is discussed. 
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1. Introduction 

A standard description of the World is usually presented in terms of the observable 
3 + 1 spatio-temporal dimensions. At the same time string theories have been developed, 
seeking to produce a consistent description of the Standard Model of physics including 
the phenomenon of gravity, which appear to be most consistent if large numbers of 
dimensions are postulated. A 26-dimensional space-time was deemed necessary for 
bosonic strings [1] and a ten-dimensional one for type-II [21 [3] and heterotic strings [1] . 
The latter theories are closely related to a mysterious theory called M-theory, which 
lives in 11 dimensions [5]. In contrast, classical physics requires 3 spatial dimensions, 
e.g. to accommodate Newton’s inverse square law as argued already by Immanuel Kant. 
Ehrenfest has shown that atoms only exhibit stable orbits in a 3-dimensional space [6]. 
These contradictions between requirements from classical and quantum physics for a 
3-dimensional space and the possibility of a theory involving higher dimensions were 
already resolved in 1926 by Klein invoking the concept of compactihcation [7]. 

In the present study the accurate results from precision measurements on molecules 
are exploited to constrain existing theories on higher dimensions. For molecular systems, 
state-of-the-art quantum level calculations of the molecular ions H^, HD+, and D^, 
all fundamental three-particle Coulomb systems, have reached the precision that the 
uncertainty becomes limited by the precision at which values of the fundamental 
mass ratios mp/rue and m^/mp are known mu [To] , although the recently improved 
determination of mp/me HI] demonstrates active progress on the experimental side. 
While experiments on the ro-vibrational spectrum of the isotopomer are still under 
way [12] the small dipole moment of the HD’*' isotopomer has enabled the accurate study 
of electric-dipole-allowed transitions in various bands [nun US]. 

In recent years great progress has also been made on the calculation of level energies 
in the neutral hydrogenic molecules. Accurate Born-Oppenheimer energies have been 
calculated for the electronic ground state of H2, HD and D2 [I6]5 ^ non- 

adiabatic interactions m HH] and relativistic and quantum electrodynamical (QED) 
corrections [IHISQ]. Now a full set of ro-vibrational level energies of all quantum states up 
to the dissociation limit is available for all three isotopomers [2T1I22] . These calculations 
on the ground electronic quantum levels were tested in experiments measuring the 
dissociation limits of H2 [23], D2 [23], and HD [25]. Further they were compared to 
experimental values for the fundamental vibrational splitting in H2 and the hydrogen 
isotopomers [26] |2Z] , to a measurement of the hrst overtone in H2 [23 [30] and D2 [22] , 
a measurement of the second overtone in H2 [^ [32] . and measurements of highly 
excited rotational levels in H2 [33] • The results from the variety of experimental 
precision measurements on both the ionic and neutral hydrogen molecules are generally 
in excellent agreement with the QED-calculations, within combined uncertainty limits 
from theory and experiment. 

The agreement between experiment and hrst-principles calculations on quantum 
level energies of molecules has inspired an interpretation of these data that goes 
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beyond molecular physics. Since weak, strong and (Newtonian) gravitational forces 
have negligible contributions to their quantum level structure, electromagnetism is the 
sole force acting between the charged particles within light molecules, and QED is the 
fully-encompassing framework to perform calculations. This makes it possible to derive 
bounds on possible fifth forces between hadrons from molecular precision experiments 
compared with QED-calculations [SH |35] . 

Theories of higher dimensions were developed with the goal to resolve the hierarchy 
problem, i.e. the vast difference of scales between that of electro-weak unification (1 
TeV) and that of the Planck scale (10^® TeV), where gravity becomes strong. By 
permitting the leakage of gravity into higher dimensions while keeping the particles 
and the three forces of the Standard Model in 3-1-1 dimensions, and invoking a 
compactihcation range for the extra dimensions exceeding 3-1-1, two different testable 
theories were phrased by Arkhani-Hamed, Dimopoulos, and Dvali [36] and by Randall 
and Sundrum [33, [38] . The mathematical formalisms of these theories can be applied to 
molecular physics test bodies, from which constraints on the compactihcation distances 
can be deduced for the former, while constraints on the brane separation or curvature 
can be derived for the latter. That is the subject of the present paper. 


2. The ADD-model 


It is the intention of the theory formulated by Arkani-Hamed, Dimopoulos, and 
Dvali [36], referred to as ADD-theory, to establish an effective Planck scale to coincide 
with the electro-weak scale by allowing gravity to propagate in extra dimensions. The 
three forces of the Standard Model, tested at very short distances in particle and atomic 
physics experiments, are considered to act locally within a 3-brane (3 spatial dimensions 
and a time dimension) embedded in a higher dimensional bulk, where gravity may act 
allowing for gravitons to escape. By this means in ADD the hierarchy problem is 
nullihed, and the so-called desert range between the electro-weak scale (Mew) of 1 TeV 
and the Planck scale (Mpi) of 10^® TeV avoided. The extension of the extra dimensions 
is necessarily limited, in the case of flat metrics considered in ADD, since experiments of 
the Cavendish-type have proven that gravity obeys the Newtonian 1/r potential beyond 
the range of 1 cm [39]. Hence, the extra dimensions are considered to be compactified 
within a range parameter While in principle the extra dimensions could exhibit 
differing range parameters, in the ADD-formalism and in the present analysis such 
difference is not made. 

The Newtonian gravitational potential may be written as: 


RN(r) = -G 


mim2 mim2 1 

r M|i r 


he 


( 1 ) 


with the Planck mass dehned as Mfi = he/G in SI units. In the following discussions, we 
adopt the natural units h = c = 1 and drop the (hc)-factor in the potentials. The extra 
n spatial dimensions proposed in the ADD theory result in a modiheation of Newtonian 
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gravity, for distances shorter than the compactihcation length range, that is consistent 
with Gauss law: 


Vadd 


77117712 1 

]\/fn+2 rf.n+1 ’ 


( 2 ) 


where the subscript 4 represents the known (3 + 1) spacetime dimensions, and 
M( 4 +„) is the full higher-dimensional Planck mass. For separations larger than the 
compactihcation length, r > the ADD potential should be in correspondence with 
the Newtonian 1/r-form 


Vadd — 


7T7i7T72 




(3) 


To be more precise, {RnY should be the compactihed volume of the extra dimensions 14, 
thus a factor of order unity might be included for a specihc compactihcation geometry. 
The Planck mass Mpi is then related to the higher-dimensional mass M( 4 +„) via: 

Ml = (4) 


Thus the fundamental mass M( 4 +„) may still be small and Mpi becomes large due to the 
compactihed volume of extra dimensions. Arkani-Hamed et al. have shown that if the 
fundamental mass is taken as Mpw one extra dimension would have a range of order 
10^° km to account for the weakness of gravity. This is incompatible with experimental 
evidence. But for two extra dimensions would be of sub-millimeter size [36] , thus at 
a range where Newtonian gravity is not hrmly tested. In our present study we will not 
set a certain energy scale, and in particular we do not assume that M 4 +„ ~ Mpw Our 
goal is to constrain Rn from molecular physics experiments without theoretical prejudice 
regarding the fundamental mass scale. 

While dealing with molecules the unit attraction of gravity can be chosen as that 
between two protons and a dimensionless gravitational coupling strength is dehned as: 


ac = Gmp/hc (5) 

Note that this particular choice of the gravitational coupling constant is equivalent to 
specifying ac = {mp/Mpi)‘^ = 5.9 x 10“^®. Then the Newtonian attraction between two 
particles consisting of A^i and N 2 protons or neutrons ( 777 ^ ~ mp is adopted) can be 
written as: 

^N(r) = -aGNiN2- (6) 

r 

From Eq. (jlD, the ADD-potential of Eq. (I2|) within the compactihcation radius r < Rn 
may be rewritten as: 

GADD(r) = -aGN,N2Rl^^, (7) 

while this potential reduces to normal Newtonian gravity 14 for the range outside the 
compactihcation length range r > 
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Figure 1. Wave functions for H 2 in the electronic ground state with z; = 0, J = 0 and 
z; = 1, J = 0, and for the HD+ z; = 0, J = 2 quantum state. 


For molecules this gravitational potential has an effect on the level energy of a 
molecular quantum state with wave function \l/(r), to be written as an expectation 
value: 


(Vadd) 


aGNiN2 


POO 

/ '^*(r)-'^(r)r‘^dr 

jRn r 

+RZ C ^-(O^^mryr^dr 


( 8 ) 


Note that the wave functions are given along a single coordinate r, i.e. the vibrational 
coordinate, that probes the gravitational forces between nucleons. Here the nuclear 
displacement is separated from electronic motion and the wave function (r) represents 
the probability that the nuclei in the molecule are at internuclear separation r. The 
hrst integral term represents the ordinary gravitational attraction, which is for protons 
8 X 10“^^ times weaker than the electrostatic repulsion, and can therefore be neglected. 
The second integral represents the effect of modihed gravity and is evaluated using 
accurate wave functions for H 2 . The wave functions of the H 2 ground electronic state 
for the n = 0 and 1 levels are shown in Fig. [H In practice, the integration is performed 
up to r = 10 A since the wave function amplitude is negligible beyond that. Also at 
shorter distances r < 0.1 A the wave function amplitude becomes negligible, for which 
reason the second integral in Eq. ([H]) converges without additional assumptions. The 
HD'*' n = 0, J = 2 ground electronic state wave function is also displayed in Fig. [1] 
showing the larger internuclear distance of the ion with respect to the neutral. 

For transitions between quantum states and T 2 , as in spectroscopic transitions 
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in molecules, a differential effect must be calculated: 


{AVAMn,Rn)) = -acNiN^K 




^n+l 




rTL+l 




(9) 


This equation represents the expectation value for a high-dimensional gravity 
contribution to transitions in molecules. Here the ADD-expectation value is written 
explicitly as a function of the two relevant parameters: the number n of extra spatial 
dimensions and the compactihcation scale Rn- From Eq. (jH]), it is clear that a stronger 
effect can be expected if the difference in wave functions of the two states Ti and T 2 is 
greater. For this reason, measurements on the dissociation limit in molecules, where Ti 
is lowest energy bound state and T 2 is the non-interacting two-atom limit at r = cxo, 
are the most sensitive probes. 

3. The Randall-Sundrum models 


Let us now consider the Randall-Sundrum scenarios, RS-I and RS-II, to approach the 
physical description of extra dimensions in an alternative manner [371 EH]. In these 
scenarios, the particles and interactions of the Standard Model are conhned in the SM- 
brane, separated by some distance yc from another (hidden) 3-brane along one extra 
dimension y. The branes and the bulk are sources of gravity that were shown to produce 
an Anti-de-Sitter metric: 

ds^ = e~‘^^^'^^rj^ydx^dx’' -\- dy'^, (10) 

where exp{—k\y\) is a so-called warp factor and k is the bulk curvature [37]. The 
warped metric differentiates the RS models from the ADD model with a flat metric 
where k = 0. Thus, the exponential warp factor in the RS scenarios solves the hierarchy 
problem alternatively, without requiring large extra dimensions as assumed in the ADD 
model. 

In the RS scenarios the modihed gravitational potential between two masses 
separated by a distance r in the SM-brane can be expressed as: 

^Rs(r) = -G^^^(l + AKs), (11) 

r 

where Ars is the correction to the Newtonian potential. Callin [ID] computed the 
potential in the framework of the RS-I scenario, obtaining for short distances: 


ARs-i(r) ~ 


4 


3 2/cyc 


Snkr 1 -|- ’ 


kr 1. 


( 12 ) 


Here one can distinguish two regimes, kyc 1 and kyc S> 1, with the result up to the 
leading order: 




ARs-i(r) ~ < 


Sirkr 


+ 


+ ... 


kyc < 1 , 

kyc > 1. 


(13) 
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It turns out that the RS potential for long distances (fcr 1) is not applicable to 
molecules and is not considered further. 

In the RS-II scenario, the hidden 3-brane is chosen to be inhnitely far {yc —?■ oo) 
from the SM-brane resulting in an effective model with a single 3-brane (SM-brane) 
in the bulk. This solution thus offers the existence of extra dimensions that do not 
require compactihcation in contrast to the ADD model. For short distances in the 
RS-II scenario, Callin and Ravndall [H] obtained 


ARs-ii(r) = ^ + ..., fcr<l, (14) 

for the RS correction. Note the correspondence of Eq. flTT)) with that of Eq. flT^ for 
kyc ^ 1, which is expected since the latter RS-I condition implies the transition to 
RS-II at inhnite brane separation. 

From these RS potential corrections, the expectation values of the leading-order 
shifts of transitions in molecules, in the short distance separation {kr 1) regime, are 
therefore: 


{AVnsik)) = acN^N^J^ 


Silk 


Ti 


Ti 






(15) 


where iF = {1 




)/(l + 


for RS-I and = 1 for RS-II. Using these 


expressions, limits on the curvature k or the brane separation yc based on molecular 
spectroscopy data can be derived. 


4. Constraints on higher dimensions from molecular data 

In the previous section, the expectation value for a higher-dimensional gravity 
contribution to a transition frequency in a molecule was presented for both ADD and RS 
approaches to higher dimensions. This expectation value is interpreted as a contribution 
to the binding energy of molecules in certain quantum states. This rationale will be 
used to derive constraints on characteristic parameters underlying the extra-dimensional 
theories, the compactihctaion range Rn for the ADD scenario and the warp factor k or 
brane separation yc for the RS scenario(s). 

In Table [T] a compilation is made of a comparison between theoretical and 
experimental values obtained in recent experiments for hydrogen neutral molecules 
and hydrogen molecular ions, and the stable isotopomers containing deuterons. Ro- 
vibrational transitions in the ground electronic state are indicated by the change 
in vibrational quantum number n, while Dq denotes the dissociation energy of the 
ground electronic state. In the Table the agreement between theory and experiment 
is represented by the combined uncertainty 5E with: 

SE = (16) 

where and ^Etheory signify uncertainties of theory and experiment. On all but 

two cases the values for 5E were found to be larger than the discrepancies between 
theory and experiment, denoted by AE = Eexp — Etheory, while the H 2 n = 0 —)■ 1 is 
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Table 1. Data from recent precision measurements of vibrational energy splittings 
as well as the dissociation energy Dq in neutral and ionic molecular hydrogen and 
their isotopomers. Adapted from Ref. [53] and updated with most recent data. AE 
represents the deviation between theory and experiment, while SE represents the 
combined uncertainties, cf. Eq. (HU). 


species 

transition 

AE (cm 

5E (cm 

Ref. 

H 2 

n = 0 —)■ 1 

0.000 24 

0.00017 

[HUT] 


v = Q^2 

0.000 4 

0.002 0 

[23130] 


71 = 0 —)■ 3 

-0.000 6 

0.002 5 

[3I1I32] 


Do 

0.000 0 

0.0012 

m 

HD 

V = 0 —)■ 1 

0.00011 

0.000 23 

[261 [27] 


Do 

0.000 9 

0.0012 

m 

D 2 

z; = 0 —)• 1 

-0.000 02 

0.00017 

[lUllT] 


z; = 0 ^ 2 

-0.000 5 

0.001 

m 


Do 

0.000 5 

0.0011 

m 

HD+ 

z; = 0 —)■ 1 

-0.000 005 2 

0.000 002 0 

[ 13112 ] 


z; = 0 —)■ 4 

0.000 009 

0.000 017 

m 


within two standard deviations {AE < 2 6E). From these results it is concluded that 
QED-theory for these molecular systems is in very good agreement with observations. 
Recent calculations by Korobov et al. |32] result in an increased discrepancy with the 
experimental results of Bressel et al. na at the level of 2.6 standard deviations, and we 
do not include the HD+ n = 0 ^ 1 values in the comparisons. 

The agreement between theory and experiment for molecular systems is now 
translated into a constraining relation for higher dimensions in the ADD framework: 


{AVABB{n,Rn)) <6E. 


(17) 


As a first example we take the measurement on the fundamental vibration in the H 2 
molecule. This is one of the most accurately measured numbers in neutral molecules, 
while also the QED-calculations for this fundamental rotationless transition are more 
accurate by an order of magnitude with respect to the absolute binding energies, because 
of cancellation of errors for non-rotating molecules [26]. Constraints on Rn can be 
derived via: 


{Rr, 


< 


6E 


(18) 


q:qNiN2A 

with A the difference in expectation values over the wave function densities between 
n = 0 and v = 1 vibrational states in the molecule: 

1 \ / 1 \ 


A = 






(19) 


The wave functions for the lowest vibrational states, in the case of H 2 and for 
J = 0, as obtained from ah initio calculations mm are plotted in Fig. [1] Since 
the wave functions are located in the same region of space the fundamental vibrational 
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R„ (m) 


Figure 2. (Color online) Limit on the compactification range i?„ as derived from the 
measurement of the fundamental vibration in the H 2 molecule |26j in comparison with 
the ADD-formalism. 


transition in the hydrogen molecule (u = 0 ^ u = 1) probes only a differential effect. 
The resulting constraints on from the measurement of the fundamental vibration in 
the H 2 molecule for the range of extra dimensions n = 1 — 8 are presented in Fig. |2l 
The sloping lines in Fig. |2] represent calculated Vadd/^^j for different n and Rn values. 
The horizontal dashed line 5E/V^ indicates limits from molecular spectroscopy. Hence, 
for certain numbers of extra dimensions n, Rn is constrained to be less than the value 
where the Vadd/Hn and SE/Vfi intersect in the graph. Constraints on Rn, obtained 
from a comparison with the fundamental vibrational transition of H 2 are presented in 
Table [2l 

The experimental as well as the theoretical results for the fundamental vibration in 
the hydrogen molecule are known to the 10“"^ cm“^ level, an order of magnitude more 
accurate than the values for the binding energies [261ET]. However, for a comparison of 
dissociation limits it is no longer a small difference along the internuclear coordinate axis 
that is probed, but the difference between the 1 A molecular scale and infinite atomic 
separation. The expectation value for the ADD-contribution to the binding energy of 
the lowest bound state in the H 2 molecule, or the Dq binding energy, is: 



( 20 ) 


By comparing to the experimental findings on Zi)o(H 2 ) [23] this leads to another set of 
constraints on Rn for n extra dimensions, which are also listed in Table El 

The method was further applied to the fundamental vibration of HD and D 2 , where 
the experimental and theoretical uncertanties are similar to those in H 2 . Although the 
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Table 2. Constraints on the size i?„ of compactified dimensions (in units of m) as 
derived from a number of molecular features: (i) the fundamental (0 —^ 1) vibration in 
H 2 , (ii) the dissociation limit Dq of H 2 , (iii) the dissociation limit of D 2 , and (iv) the 
(4-0) R(2) ro-vibrational transition in HD+. The constraints are derived within the 
ADD-framework assuming that n extra dimensions are of equal size. The corresponding 
higher-dimensional Planck length i?pi,( 4 +„) (in units of m) and Planck mass 
(in units of GeV) are also tabulated, where the smallest values for i?pi^( 4 +„) and the 
highest value for M( 4 +„) is taken from the examples. 


n 



Rn 


Rpi,(4-|-n) 

M(A+n) 


H2 (1-0) 

H2 Do 

D2 Do 

HD+ (4-0) 

(m) 

(GeV) 

2 

2.2 X 10^ 

1.0 X 10^ 

4.8 X 103 

2.8 X 103 

2.1 X 10-16 

9.3 X 10-1 

3 

7.7 X 10-1 

1.9 X 10-1 

1.2 X 10-1 

1.0 X 10-1 

3.0 X 10-13 

6.5 X 10-2 

4 

1.1 X 10-3 

8.5 X 10-^ 

5.9 X 10-4 

7.0 X 10-4 

1.8 X 10-14 

1.1 X 10-2 

5 

3.3 X 10-3 

3.2 X 10-3 

2.4 X 10-3 

3.1 X 10-3 

5.8 X 10-14 

3.4 X 10-3 

6 

3.4 X 10-6 

3.7 X 10-6 

2.9 X 10-6 

3.0 X 10-6 

1.4 X 10-13 

1.4 X 10-3 

7 

6.9 X 10-1' 

7.8 X lO-i" 

6.4 X lO-"^ 

6.3 X 10-1' 

2.8 X 10-13 

7.1 X 10-4 


heavier masses of the isotopomers improve the constraints obtained from H 2 , as expected 
from Eq. IfTSll . the HD and D 2 fundamental vibration constraints are still less stringent 
compared to that from the H 2 dissociation limit. The results obtained for D 2 dissociation 
energy [21] lead to the tightest constraints on from the neutrals as listed in Table [21 
which scale by a factor ( 1 / 4 )^/"' relative to H 2 due to the mass difference. 

The experimental accuracy for the molecular ion transitions is an order 

magnitude better than the corresponding neutral molecule system that stems mostly 
from the possibility of trapping the ionic species. The theoretical calculation for the 
three-body HD’*' level energies is also more accurate than those of the neutral molecular 
hydrogen. However, the internuclear separation of HD”*^ (~ 1.1 A) is greater than that 
of neutral hydrogen molecules (~ 0.76 A) as shown in Fig. [H Thus the neutrals are 
inherently more sensitive as the wave functions probe shorter internuclear distances 
compared to their ionic counterparts. The constraints for derived from the HD’*' 
(n = 0, J = 2 —)■ n = 4, J = 3) ro-vibrational transition from Koelemeij et ah [I3] are 
listed in Table [21 In the table, the constraints from D 2 Dq are the most stringent for 
n = 4, 5, 6 extra dimensions while the constraints from HD+ are the most constraining 
for n = 2, 3, 7. The higher-dimensional Planck mass and corresponding Planck 

length i?pi,( 4 +n) derived from the tightest Rn constraints obtained in this study are also 
listed in Table [21 

Similarly, we derive constraints pertaining to corrections in the RS scenario with 
one extra dimension, and the combined uncertainty 6E for a specihc molecular transition 

(AVnsik)) < 6E. (21) 

Using the combined uncertainties for the Do(D 2 ) study, we present constraints for the 
RS schemes. For the RS-I scenario in the short distance {kr 1) regime, we obtain 
constraints for the brane separation of |/c < 1 x 10^® m in the limit kyc 1 using 
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Eq. (fT^ . In the limit kyc ^ 1 in Eq. fll3p . a constraint for the inverse of the curvatnre 
of l//c < 2 X 10^® m is obtained. For the RS-II model, we obtain constraints for the 
inverse of the curvatnre 1/fc < 2 x 10^® m for fcr -C 1 from Eq. (1T5|) . 

5. Comparison with other constraints 

The constraints obtained from molecular systems probe length scales in the order of 
Angstroms. This complements bounds probing subatomic to astronomical length scales 
obtained from other studies using distinct methodologies. Length scales of several 
hundred nanometers to microns are probed in Casimir-force studies using cantilevers [13] 
or atomic-force microscopy [S]. The micrometer to millimeter range are probed in 
torsion-balance type experiments, with the tightest constraint obtained by Kapner et 
al. j45] for a single extra dimension of Ri < 4.4 x 10“^ m. The centimeter to meter 
separations are accessed by Cavendish- or Eotvos-type investigations in the laboratory, 
while astronomical scales can be probed in satellite or planetary orbits that also serve 
to constrain the universality of free fall and deviations from the gravitational inverse- 
square law [39]. Constraints for the RS-theories are obtained by lorio [16] using data 
from orbital motions of satellites or astronomical objects, with the tightest constraint 
for the inverse of the curvature of 1/fc < 5 m obtained from the motion of the GRACE 
satellite. The latter constraint is in the fcr S> 1 regime of Eq. ffT4|) and probes a different 
distance range to that of molecules {kr <C 1). 

Precision spectroscopies of hydrogen [471148] and muonic atoms [49l |50] have been 
interpreted along the same lines in terms of the ADD-model ISB resulting in typical 
constraints of < 10“® m. The interpretation is not straightforward because of 
the proton size puzzle [52]; in fact, the argument has been turned around, where the 
existence of extra dimensions are instead invoked as a possible solution to the puzzle [53] . 
In the treatment of atoms, some assumptions had to be made on the wave function 
density at r = 0, typical for the s-states involved, causing problems in calculating 
the second integral of Eq. ([8]) over the electronic wave function that has a significant 
wave function amplitude at r = 0 in atoms. Note that these difficulties are absent in 
molecules, as the molecular wave function probes the 0.1-5 Angstrom distance range. 

To probe length scales in the subatomic range, one is ultimately limited by the 
increasing contributions from nuclear structure and the strong interaction, e.g. in 
neutron scattering studies [54] • In contrast to QED calculations, the most accurate 
lattice-QCD calculation of light hadron masses only achieves relative accuracies in the 
order of a few percent. Nevertheless, the smaller nucleon size presents higher sensitivity 
to effects of ADD-type interactions, and constraints for the size of extra dimensions 
may be extracted. The general method for molecules presented here may be applied 
to a comparison of ab initio lattice-QCD calculations with the measurements of light 
hadron masses. The corresponding QCD test probes length scales of the size of a nucleon 
at ~ 10“^^ m. The ab initio calculations of Diirr et al. [55] for the nucleon mass are 
estimated to be accurate to around 50 MeV/c^ while the experimental mass values are 
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accurate to 20 eV/c^. The calculated nucleon mass, with 3% relative accuracy, is the 
isospin average of proton mp and neutron masses, while is known experimentally 
to be 0.1% smaller than Constraints based on experimental nucleon masses and 
QCD theory have not been explored, but we produce here a rough hrst estimate by 
assuming that the three constituent quarks each have an effective mass that is ^/s of 
the nucleon mass, and have separation distances ~ r^. Analogous to Eq. ffT7|) for QED 
interactions in molecules, the expectation value for an ADD contribution on the mass 
of the proton can be written as Vadd(p)/c^ < 5mp, yielding a bound for the case of 7 
extra dimensions of i ?7 < 2.4 x 10 m. 

In high-energy particle collisions, higher-dimensional gravitons may be produced 
that could escape into the bulk, leading to events with missing energy in (3-1-1)- 
dimensional spacetime [36l [56| IFf] . Based on this premise of an energy loss mechanism 
the phenomenology of the SN 1987A supernova was investigated, imposing limits on 
extra dimensions of i ?2 < 3 x 10“® m, i ?3 < 4 x 10“^ m, and < 2 x 10“® m [58] . 
Similarly, from a missing energy analysis of proton colliding events at LHC, a constraint 
for i ?2 < 3.2 X 10“^ m can be extracted from the = 1-93 TeV bound for n = 2 given 

in Ref. [59] . For comparison, the Planck energy scale in (4 -|- n) dimensions in Table [2] 
turns out to be in the range between 1 — 1000 MeV, but are derived from a completely 
independent methodology. Also for n > 2, the bounds derived from LHC are nominally 
more stringent than those from molecules. However, additional assumptions beyond the 
ADD potential in Eq. ((2|), e.g. the fundamental quantization of gravity, the existence 
and propagation of gravitons in (4 -|- n) dimensions and postulating the existence of 
massive new particles, are necessary for an effective theory [561 EZj to interpret the 
LHC missing energy signals. Such assumptions are not needed for the molecular physics 
bounds, which are not sensitive to physics at very short distances. 

6. Conclusion and outlook 

The alternative approaches to constraining compactihcation radii for extra dimensions, 
partially surveyed here, are all complementary as they probe different length and 
energy scales. Some approaches serve to produce tighter limits, however, often at the 
expense of additional assumptions. In the present study, a constraint is derived on 
compactihcation scales of extra dimensions from precision measurements on molecules, 
leading to straightforward interpretations. Molecules, in particular the lightest ones 
as neutral and ionic molecular hydrogen, exhibit wave functions representing the 
internuclear distances, with amplitudes conhned to the range 0.1 — 5 A. Current state- 
of-the-art experiments on neutral molecular hydrogen determine vibrational splittings 
on the order of 10“^ cm“^, or 3 MHz [26] . Since the lifetimes of ro-vibrational quantum 
states in H 2 are of the order of 10® s [60], measurements of vibrational splittings on 
the order of 10^^ Hz could in principle be possible at more than 20-digit precision, 
which leaves room for improvement “at the bottom” of over 10 orders of magnitude, 
if experimental techniques can be developed accordingly. Similar improvements in 


Constraints on extra dimensions from precision molecular spectroscopy 


13 


theory would make these molecular systems an ideal test ground for constraining or 
detecting higher dimensions, as well as fifth forces [3l]. After having performed a 15- 
digit accuracy calculation on Born-Oppenheimer energies [16], calculations of strongly 
improved accuracy have just been published [6T|, while improved calculations of non- 
adiabatic corrections are underway [62] • Immediate improvements, based on existing 
technologies, on the experimental accuracies of the dissociation limits in the neutral 
hydrogen and its isotopomers [63] and the spectroscopy of HD+ [121 Ell ES], were 
discussed recently. 
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